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We present a search for hypothetical vector-like quarks in pp collisions at yfs = 1.96 TeV. The 
data were collected by the DO detector at the Fermilab Tevatron Collider and correspond to an 
integrated luminosity of 5.4 fb" 1 . We select events with a final state composed of a W or Z boson 
and a jet consistent with a heavy object decay. We observe no significant excess in comparison to 
the background prediction and set limits on production cross sections for vector-like quarks decaying 
to VF+jet and Z+jet. These are the most stringent limits to date for electroweak single vector-like 
quark production at hadron colliders. 

PACS numbers: 12.60.-i, 13.85.Rm, 14.65. Jk 



The standard model (SM) of particle physics, despite 
its many successes in accurately describing interactions 
below the TeV scale, is known to suffer shortcomings 
at higher energy scales. A wide range of theories have 
therefore been proposed to describe phenomena at the 
TeV scale and beyond, among them warped extra di- 
mensions universal extra dimensions Q, and little 
Higgs j| models. There exist particular realizations 
of each of these theories that predict the existence of 
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vector- like quarks [4-6], massive particles which share 
many characteristics with SM quarks. The notable ex- 
ception is that the right-handed and left-handed compo- 
nents of vector-like quarks transform in the same way 
under SU(3) x SU{2) x U(l), hence their name. 

Previous searches for strong pair production of vector- 
like quarks at the Fermilab Tevatron Collider have ex- 
cluded masses of up to 338 GeV at the 95% C.L. 0. 
However, recent years have seen the development of the- 
oretical scenarios in which corrections to SM quark cou- 
plings arising from their mixing with vector-like quarks 
cancel out, permitting the two types of quarks to mix 
to a degree unconstrained by precision electroweak mea- 
surements and 6-factory results. Such large mixing with 
no introduction of anomalous couplings allows for single 
production of vector-like quarks at hadron colliders via 
the weak interaction 0] ■ Diagrams for single electroweak 
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FIG. 1: s-channel (a) and f-channel (b) electroweak produc- 
tion of a vector-like quark Q at the Tevatron. 

vector-like quark production are shown in Figure [TJ 

Generally, electroweak couplings to SM quarks are set 
by the parameter k q q: 

v 

K qQ = K q Q, (1) 

m Q 

where v is the vacuum expectation value of the SM Higgs 
field, mg is the mass of the vector-like quark, and R q Q 
is the coupling strength. Provided this coupling is not 
too small, a search for singly produced vector-like quarks 
can benefit from the lower kinematic threshold compared 
to pair production and stands to improve the mass limit 
considerably. 

In this Letter, we present a search for singly produced 
vector-like quarks using data corresponding to an inte- 
grated luminosity of 5.4 fb" 1 collected by the DO detec- 
tor at the Fermilab Tevatron Collider. We consider a 
model [H in which there are two doublets of vector-like 
quarks that couple to the first generation of SM quarks. 
Vector-like quark final states are characterized by either a 
W or Z boson and at least two jets, one of which results 
from the decay of the vector-like quark, while another 
is produced in association with the vector-like quark at 
the primary vertex. We select events in which the vector 
boson decays to either electron(s) or muonfs). 

The DO detector is described in detail in [8| . The region 
of the DO detector closest to the interaction point con- 
tains a central-tracking system consisting of a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT), 
both of which are located within a 2 T superconducting 
solenoidal magnet. Hits in these two detectors are used to 
reconstruct tracks from charged particles. Surrounding 
the two tracking subdetectors are liquid-argon and ura- 
nium calorimeters, both electromagnetic and hadronic, 
which have a central section (CC) covering pseudorapidi- 
ties Q \ri\ < 1.1, and two end calorimeters (EC) that 
extend coverage to \r/\ f=s 4.2, with all three housed in 
separate cryostats. Electrons are identified as isolated 
energy clusters in the electromagnetic calorimeter, with 
a shape consistent with that of an electron, matched to 
a track in the inner detector. Jets are reconstructed in 
the calorimeters using the iterative midpoint cone algo- 



rithm [10] with a cone of radius 1Z — 0.5 in rj — <p space. 
An outer muon system, providing coverage for \rj\ < 2, 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T toroids, followed by two 
similar layers after the toroids. Muons are identified via 
a combination of reconstructed tracks in both the outer 
muon system and the central tracking system, and are re- 
quired to be isolated from both tracker and calorimeter 
activity. 

Vector-like quark signals and background processes 
which include electrons or muons from W or Z boson 
decays are modeled using Monte Carlo (MC) simula- 
tion. Signal samples are generated using MADGRAPH 
with CTEQ6L1 [I3|_parton distribution functions, LO 
cross sections from |4( and vector-like quark resonance 
widths calculated with bridge [l3|. Subse quen t parton 
shower evolution is generated with PYTHIA [14|. We set 
k„d = 1, K u u = \[2 and kdu = K-dD = 0, such that down- 
type vector-like quarks decay exclusively to Wq and up- 
type quarks decay exclusively to Zq. This choice, while 
affecting the production rates and therefore the excluded 
vector-like quark masses, does not affect the cross sec- 
tion limits themselves. We generate up- and down-type 
vector-like quark signals with masses between 280 GeV 
and 700 GeV. Backgrounds from tt, W — > iu+jets, and 
Z — > tt+jets production are modeled using ALPGEN [l5j |. 
also interfaced to pythia. For these events, we use a 
matching procedure to avoid double counting partons 
produced by ALPGEN and those added by the shower- 
ing in pythia (lij . Diboson samples are generated with 
pythia, and single top quark production is modeled us- 
ing the COMPHEP [13] generator. All MC samples are 
passed through a GEANT [l8[ simulation of the DO detec- 
tor and are overlaid with data events from randomly cho- 
sen beam crossings to simulate the effect of multiple pp 
interactions and detector noise. MC events are then re- 
constructed using the same software that is used for data 
reconstruction. The VF+jets and Z+jets samples are 
normalized to the leading-log cross section reported by 
ALPGEN times a fc-factor calculated by MCFM [l9j]. The 
tt samples are normalized to a next-to-next-to-leading or- 
der (NNLO) cross section for m t = 172.5 GeV [13], and 
the diboson samples (WW, WZ and ZZ) to the NLO 
cross section predicted by MCFM pi) . All simulated 
events are corrected for differences in trigger and recon- 
struction efficiencies between data and simulation. 

We conduct the search in two channels, correspond- 
ing to vector-like quark decays to {W — > lv)+jet and 
to (Z — s- ££)+}et. In the IV+jets channel, we initially 
select events which have passed a single lepton trigger 
and contain exactly one electron (\rj\ < 1.1) or muon 
(I77I < 2.0) with transverse momentum pr > 20 GeV, at 
least two jets with px > 20 GeV and \ij\ < 2.5, and miss- 
ing transverse energy corrected for the momentum 
of any muons in the event, greater than 15 GeV. We also 
require 2Mjf + $ T > 80 GeV, where is the trans- 
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verse mass of the lepton-^ T system, in order to suppress 
the multijet background. In the Z+jets channel, we ini- 
tially select events with exactly two electrons or muons 
in addition to at least two jets, all with same pt thresh- 
olds as used in the VF+jets channel. Events are selected 
by a mixture of single lepton, dilepton and lepton plus 
jets triggers. Due to the low rate of background without 
real Z bosons, electrons are also accepted in the endcap 
electromagnetic calorimeters (1.5 < \rj\ < 2.5). The two 
lcptons are required to have an invariant mass between 
70 and 110 GeV, i.e. consistent with that of a Z boson. 
Additionally, we require lfi T < 50 GeV, as this channel 
contains only instrumental sources of missing transverse 
energy. 

The largest physics background to the single lepton 
channel is W(— > £v)+jets production, with smaller con- 
tributions from Z+jets, tt, diboson, and single top quark 
processes. The main instrumental background arises 
from multijet events in which one of the jets is misiden- 
tified as a high-pr isolated lepton in the detector. We 
model this background using data events which fail the 
calorimeter shower shape requirements for the electron 
or the isolation requirements for the muon, but pass all 
other selection criteria. In the single electron channel, 
we estimate the relative fraction of real electrons from 
W boson decays and misidentified electrons from jets by 
determining the efficiencies for each type of event to pass 
a tighter selection. These efficiencies are calculated in 
Z — > ee data for real electrons and in a ~$ T < 20 GeV 
sample for misidentified electrons. In the single muon 
channel, we scale the events failing the muon isolation 
requirement to match the number of events obtained af- 
ter subtracting the expected number of real W bosons 
satisfying the event selection from the number of events 
observed in data. 

In the dilepton channel, (Z — > ££)+jets events domi- 
nate the SM background. In order to correct the Monte 
Carlo for the small trigger inefficiency in the data, we 
apply a global normalization determined using a pure 
Z — >• l£ sample. The ratio between the data yield in 
this sample and the predicted yield from the inclusive 
Z boson cross section multiplied by the branching ra- 
tio is associated with the overall trigger efficiency, and 
the Monte Carlo distributions are scaled by this ratio. 
The multijet background is modeled using data events in 
which both leptons fail quality criteria (for electrons) or 
isolation criteria (for muons), and is normalized to the 
difference between data and MC in the dilepton mass 
window between 40 and 70 GeV. 

The signal contains events with a W or Z boson and 
jets, all with high transverse momentum, as the de- 
cay products of a high mass resonance. We apply sev- 
eral kinematic selection criteria to select events of this 
type and to minimize the contributions from SM back- 
ground processes. In the single lepton channel, we re- 
quire lepton pt > 50 GeV, highest jet pr > 100 GeV, 



TABLE I: Predicted number of background events including 
total uncertainties and observed number of data events after 
final selection. 



Process 


Single lepton sample 


Dilepton sample 


Multijet 


47.7 ±4.7 


< 0.1 


Z+jets 


39.9 ±7.4 


262 ±45 


W+jets 


901 ± 159 


0.3±0.2 


Top 


193 ±24 


0.57 ±0.06 


Diboson 


38.6 ±3.8 


8.3±0.7 


Background sum 


1220 ±161 


271 ±45 


Data 


1175 


285 



ijjrp > 40(50) GeV for the muon (electron) channel, and 
that the separation in azimuthal angle <fi between the lep- 
ton and the ~$ T be less than 2. As the signal contains 
real W bosons, we additionally require Mjf < 150 GeV. 
Finally, we exploit the relationship between the lepton 
charge and the rj of the jet with the second-highest px in 
the signal topology. This jet, which we assign to the SM 
quark produced in association with the vector-like quark, 
is emitted, in general, into one of the forward regions of 
the detector, and its direction is strongly correlated with 
the charge of the produced vector-like quark, and thus 
also with the charge of the lepton from its decay. We 
therefore require Qg x iy 2 > 0, where Qi is the lepton 
charge and r/j 2 is the rj of the jet with the second-highest 
Pt in the event. This selection is efficient for the signal 
(« 85%), while reducing the SM background by roughly 
a factor of two. 

In the dilepton channel, we similarly select events with 
properties characteristic of a heavy resonance decay to 
a Z boson and a jet. We require the px of the dilepton 
system to be greater than 100 GeV, the spatial separation 
of the two leptons in rj — <f> space (A7£) to be less than 
2.0, and the leading jet px to be greater than 100 GeV. 

Table U displays the observed number of data events, 
along with estimated background yields, after the final 
event selection in the single lepton and dilepton chan- 
nels. Fig. [2] shows the reconstructed vector-like quark 
transverse mass, defined as 

(M^) 2 = ^p\ w +M& + PTjl ) 2 - {p Tw +P T]1 )\ (2) 

for the single lepton channel, where Mw = 80.3 GeV is 
the mass of the W boson and pr„ refers to the transverse 
momentum of the leading jet in the event. Also shown in 
Fig. [2] is the vector-like quark mass in the dilepton chan- 
nel, reconstructed as the invariant mass of the dilepton 
± leading jet system. 

Major sources of systematic uncertainty common to 
both analyses include modeling of W-7Z±jets back- 
grounds (15%); cross sections for tt (10%), diboson (6%), 
and W/Z+jets (6%) production; and jet energy scale 
and resolution (l%-5%). Major systematic uncertain- 
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FIG. 2: (a) Vector-like quark transverse mass and (b) vector- 
like quark mass for the single lepton and dilepton channels, 
respectively. Distributions for signal processes Qd — >• Wq 
and Qu — > Zq are normalized to the integrated luminosity of 
the data and include detector acceptance and reconstruction 
efficiencies. These assume k„d = 1 and k u u — y2. 



ties unique to the single lepton analysis include inte- 
grated luminosity (6.1%) [22(, lepton identification effi- 
ciencies (3%), high-pT muon modeling, and trigger mod- 
eling (1%). A systematic uncertainty of 5% is assigned to 
the global background normalization applied to the dilep- 
ton MC samples. Systematic uncertainties on multijet 
normalization are 6.5%-T00%, depending on the chan- 
nel. These do not have a large effect on the overall back- 
ground prediction, as the estimated multijet background 
is small after the final selection. 

We observe no significant excess of data over the back- 
ground prediction in either channel, therefore we set lim- 
its on vector-like quark production cross sections. We 
employ a modified frequentist approach using the DO like- 
lihood fitter 23[ that incorporates a log-likelihood ratio 
(LLR) test statistic j24[. We calculate confidence levels 
for the signal plus background (CL s+ h) and background- 
only (CLb) hypotheses by integrating the (LLR) dis- 
tributions obtained by generating pseudo-experiments 
using Poisson statistics. Systematic uncertainties are 
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FIG. 3: Excluded production cross section for a vector-like 
quark Q decaying to Wq as a function of mq, compared to 
LO predictions of vector-like quark production with different 
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FIG. 4: Excluded production cross section for a vector-like 
quark Q decaying to Zq as a function of mQ, compared to 
LO predictions of vector-like quark production with different 



treated as uncertainties on the expected number of sig- 
nal and background events. The 95% confidence level 
limit is defined as the signal cross section for which 
CL S = CLg+b/CLi, is 0.05. Using as discriminant vari- 
ables the vector-like quark transverse mass for decays to 
Wq and the vector-like quark mass for decays to Zq, we 
obtain 95% C.L. limits on vector-like quark production 
cross sections. These are shown in Figs. [3] and HI along 
with leading order theoretical predictions for two differ- 
ent scenarios. For the case of k u d = 1 and k u u = \/2 
with no coupling to the down quark, we exclude masses 
below 693 GeV for a vector-like quark decaying exclu- 
sively to Wq and masses below 551 GeV for a vector- like 
quark decaying exclusively to Zq at the 95% C.L. For 
an alternate scenario defined by Rdu = 1 and kdD = V% 
with no coupling to the up quark, the corresponding mass 
limits are 403 GeV and 430 GeV, respectively. 
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In summary, we have presented a search for single 
vector-like quark production at the Tevatron in the 
TK+jets and Z+jets final states. The observed data are 
consistent with the background expectation and limits 
on vector-like quark cross sections are derived. We ex- 
clude vector-like quark masses below 693 GeV for decays 
to Wq and masses below 551 GeV for decays to Zq at 
95% C.L. assuming vector-like quark - SM quark mixing 
parameters k u d — 1 and k u u = V% along with 100% 
branching ratios to the respective final states. 
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